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The problems of nonlinearity and high dimension have so far prevented a complete solution of the 
control of turbulent flow. Addressing the problem of nonlinearity, we propose a flow control strategy 
which ensures that the energy of any perturbation to the target profile decays monotonically. The 
controller's estimate of the flow state is similarly guaranteed to converge to the true value. We 
present a one-time off-line synthesis procedure, which generalises to accommodate more restrictive 
actuation and sensing arrangements, with conditions for existence for the controller given in this 
case. The control is tested in turbulent channel flow (Re T = 100) using full-domain sensing and 
actuation on the wall-normal velocity. Concentrated at the point of maximum inflection in the mean 
profile, the control directly counters the supply of turbulence energy arising from the interaction of 
the wall-normal perturbations with the flow shear. It is found that the control is only required for the 
larger-scale motions, specifically those above the scale of the mean streak spacing. Minimal control 
effort is required once laminar flow is achieved. The response of the near-wall flow is examined in 
detail, with particular emphasis on the pressure and wall-normal velocity fields, in the context of 
Landahl's theory of sheared turbulence. 



I. INTRODUCTION 

The return of turbulent wall flow to the laminar state is a problem with diverse and important applications, 
such as those in the aeronautics, shipping and oil industries. A comprehensive solution of the flow control 
problem still faces serious challenges. The main obstacles to a complete flow control theory are that the 
governing equations are nonlinear and of infinite dimension. This paper aims to address the problem of the 
nonlinearity, and in doing so, improve the understanding of the physical processes. As a consequence, some 
possible approaches to the problem of high dimensionality present themselves. 

In addition to the nonlinearity, flow control strategies must deal with model uncertainty and the exogenous 
disturbances such as vibration and free-stream disturbances that arise in realistic applications. This would 
suggest the use of closed-loop control strategies, which involve the feeding back of measurements of the 
system output into current and future control decisions, over open-loop strategies, where the actual system 
output is not compared to the desired output. Of the two classes, closed-loop control, offers superior 
robustness characteristics in the face of modelling error, state uncertainty and exogenous disturbance. 

Hoc control theory has been notably successful in providing good control performance for systems with 
the large class of bounded uncertainties, nonlinearities and exogenous disturbances^^. Importantly for us, 
the Tioo theory can be generalised to cope with the nonlinearity we face with the Navier-Stokes equations. 

Modern control methods typically make assumptions about the kind of model error or disturbance that 
is present in the system. H2 or optimal control assumes Gaussian state and measurement disturbances and 
Hoo control typically gives stability guarantees for model errors up to a certain bound. The Navier-Stokes 



a ) a.s.sharma @ shcfficld.ac.uk; http://www.sheffield.ac.uk/acse/stafE/sharma 

b ) j.morrison@imperial.ac.uk 



2 



nonlinearity is neither stochastic nor bounded. However in a closed or periodic domain, it is well known to be 
conservative with respect to the perturbation energy. In the following, we exploit this fact using the passivity 
theorem. The resulting control gives global stability guarantees (and consequently relaminarisation) for 
the discretised, controlled Navier-Stokes flow, where actuation and measurement requirements are met. 
This approach has been proposed 8 -^ and outlined^ in previous work. Results of our approach applied to 
simulations of turbulent flows are presented here for the first time, with application to turbulent channel 
flow. 

From a control perspective, this work advances on previous work in at least three important respects. 
First, the current approach is demonstrated for a turbulent, three-dimensional flow. Second, it offers a 
constructive synthesis procedure. Third, it provides limits on the turbulent energy production where global 
stabilisation is not possible due to insufficient actuation or sensing. 

With reference to known existing theory, the work also describes the physics behind the controller action 
and explains why a linear control strategy is always sufficient to attenuate turbulence. 



A. Modern flow control 

A comprehensive review of modern feedback flow control is available in the paper by Kim and Bewley^. 
Further physical insight is given in the recent paper by Kim«. %2 ( "optimal" ) and %oo ( "robust" ) designs 
have been applied to the linearised transition delay problem for particular wavenumber pairs by Bewley 
and Liu 85 . The optimal control approach was tried first by Joslin 86 . Joshi et al. SJ have proposed a simple 
controller design using classical methods which was used with some success in the stabilisation of infinitesimal 
and finite-amplitude disturbances. Hogberg et al£^ have demonstrated that linear feedback control can be 
used to increase the threshold perturbation amplitudes for transition to occur. In another work, Hogberg 
et alS^ presented a gain-scheduling approach which relaminarised turbulent flows in all instances tested, 
however no proof of global stability was offered. Further results on linear flow control methods are available 
in the book by Amo and Krstic 90 . The paper by Fukagata et al£^ investigates interior forcing targeting the 
Reynolds stress terms directly. 

These H2 and Hoc linear flow control strategies have been designed to delay transition, by preventing 
the flow from leaving the regime of small perturbations to the desired laminar flow. Given the likelihood 
in practical situations of large excursions due to transients, exogenous disturbances, or model error, the 
assumption of small perturbations seems inappropriate. 

In contrast, the remarkable feature of our approach is that a linear synthesis problem provides a global, 
nonlinear stability result. An implicit consequence of the stability result is the convergence of the estimation 
problem. For the first time, the controller's internal representation of the flow is guaranteed to converge to 
the true value. In contrast to the 1-Li setting, the control and estimation problems are inexorably coupled 
and the "separation principle" 7 - does not apply. 

Perhaps most pertinent to the current approach is the work described by Balogh, Liu and Krstio^. Their 
approach uses a Lyapunov stability argument (comparable to the passivity argument used here) to prove 
the existence of a globally stabilising linear control, using relatively realistic tangential boundary actuation 
and shear stress measurements. A proof of global asymptotic stability is offered for the two-dimensional 
case and it is claimed that the method will work in three dimensional flows. Their proof is applicable only 
at "sufficiently low" Reynolds numbers. 



B. An overview of the passivity approach 

Conceptually, the method of control that we propose is simple, although the details of the synthesis are 
more complicated. The idea is to eliminate perturbations to a laminar flow, which we take as the operating 
point of the system. We do not model fluctuations about the mean profile. As such, the difference between 
the turbulent mean and the laminar profile would be included in the perturbations that we wish to decay. 

The flow at any particular wavenumber experiences a coupling from other wavenumbers, via the nonlinear 
convective term, that we consider as a forcing. This nonlinear term does not produce or destroy perturbation 
energy. Consequently, if by means of feedback control, the system at every wavenumber can be made to 
dissipate energy, then the system as a whole will still dissipate energy, even with the nonlinearity. This 
objective involves only the linear system at each wavenumber and can therefore be solved using linear 
synthesis methods. 
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There is a direct analogy from circuit theory. The energy in any circuit made up of passive components 
like capacitors, resistors, inductors, will always decay, in the absence of a non-passive component such as a 
battery. 

In the case of turbulent flow, relaminarisation could be achieved simply by introducing a large dissipative 
term, analogous to a resistor. Such an approach of adding more viscosity cannot generalise to other types 
of actuation, and so is not followed. 

Hcnningson & Reddy 93 showed that non-normality of the system matrix governing perturbations to the 
laminar flow solution is a necessary condition for subcritical turbulence; hence, imposing linear stability and 
destroying the non- normality by means of feedback control provides a sufficient condition for laminar flow. 

That this condition is linear greatly simplifies the synthesis procedure, allowing use of the superposition 
principle, so that we can find the controller at each Fourier mode separately, without sacrificing the nonlinear 
stability result. The controller may be found once, off-line, since it applies at all states at a given Reynolds 
number. Since the controller is linear, the most complicated mathematical operation it must perform on- 
line is a small number of matrix multiplications at each time step. This stands in contrast to nonlinear 
adjoint-based approaches and has obvious beneficial implications for eventual practical implementation. 
With the further application of model- or controller-reduction methods, the computations at this step could 
be reduced further. 

Although many results for finite-dimensional systems in the control literature have analogues in the 
infinite-dimensional setting, the finite-dimensional theory is usually much simpler and is computationally 
tractable. This is particularly true for the case of the Navier-Stokes equations, where proofs of even rudi- 
mentary properties remain elusive^. When the system equations are discretised, as is typically done for 
practical control problems, the finite-dimensional control theory also provides a controller synthesis proce- 
dure. The procedure we use is described in the literature^ and involves a number of transformations of the 
linearised equations, then the solution of two Algebraic Riccati Equations (AREs). Whether solutions to 
these AREs exist depends on the actuation and measurement capabilities available to the control algorithm. 
Because the conditions for existence of AREs are well understood and easily checked, this problem formu- 
lation can inform the system designer about the suitability of proposed sensor and actuation arrangements 
at the design stage. 

The passivity requirement is a conservative control strategy and may not be achievable with restricted 
actuation or sensing capability. To handle this case, the synthesis methodology presented iteratively ap- 
proaches the ideal case and offers bounds on the worst-case perturbation energy production permitted by 
the controlled flow. Because we may relax the control objective of passivity, we envisage that the current 
framework may be applied to more restrictive types of body forcing. 

For this study, full-field volume forcing and measurement of wall-normal velocity only is applied to periodic, 
turbulent channel flow at Re T = 100. This simplified forcing and measurement arrangement is chosen to 
avoid confusion between the relative importance of various physical effects and the choice of any particular 
type of actuation. 



C. Linear processes and structure in wall turbulence 

The importance of linear mechanisms in turbulence has been understood since Batchelor and Proudman 96 
put forward their theory of rapid distortion (for a review see Hunt and Carruthers2£). Phillips 9 ^ has suggested 
that, in shear flows, the Reynolds stresses arise from the direct interaction between the turbulence and the 
mean shear rather than a result of indirect, nonlinear interactions. This implies that linear control schemes 
taking advantage of these mechanisms may be successful in attenuating turbulence. 

Lee, Kim and Moin 9 - have shown how many of the important dynamical processes are captured by rapid 
distortion theory (RDT) and comparison with direct numerical simulation (DNS) of sheared homogeneous 
turbulence shows that, despite its linear approximations, RDT retains "the essential mechanism for the 
development of turbulence structures in the presence of high shear rate typical of the near-wall region in a 
turbulent shear flow" . The qualitatively similar theories of Ellingsen and Palm 100 and Landah l 101 i 102 offer 
an explanation for the importance of three-dimensional disturbances in inviscid, parallel shear flow: the 
latter offers an explanation for the formation and lift-up of near-wall streaks by an algebraic instability and 
the present work is interpreted in the context of Landahl's ideas. The importance of the theory comes from 
the fact that it is both linear and is based on an analysis of the inviscid Rayleigh equation, suggesting that 
it is relevant to dynamical processes at any Reynolds number. 

Non- normality appears as a recurrent theme in the literature on transition (Schmid and Henningson 103 , 
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Schmidt). Butler and FarrelliS 5 . investigated initial conditions which are capable of the greatest energy 
growth (which they call "global optimal perturbations") at a given Reynolds number. They found that, 
in plane Poiseuille flow, the global optimal perturbation consists of a pair of counter-rotating streamwise 
vortices, though other modes should not be underestimated: in particular, oblique modes grow less, but 
faster. It remains an open question how important the initial condition problem per se is in flows that are 
already turbulent. 

Farrell and Ioanno u 106 ' 107 have suggested that the linearised Navier-Stokes equations in plane channel 
flow under stochastic forcing can exhibit behaviour reminiscent of the streamwise vortices and streaks char- 
acteristic of turbulent flow. Kim and Limi2£ demonstrated in simulations of turbulent channel flow that 
the turbulence decays without the term coupling the wall-normal vorticity and the wall-normal velocity 
in the linearised Navier-Stokes equations. Henningson and Reddy^ have shown that non- normality is a 
necessary condition for sub-critical transition, i.e. the linearised Navier-Stokes equations must have either 
exponentially growing modes or transiently growing solutions for transition to occur. 

The prevalence of streaks and quasi-streamwise vortices in near- wall turbulence has been known for some 
time (Kline et a/j^, Kim et alM&, Robinsor>iii) , though which feature causes the other is still a subject 
of discussion (see for example Chernyshenko^) . More recently, Kim & Adrian^, Ganapathisubramani et 
alJ^, Hutchins & Marusio^ and Guala et al*^ have shown the importance of very large scale motions 
(VLSMs or "superstructures") which carry approximately half the Reynolds stress. Recently, they have 
been shown to appear on very rough surfaces also (Birch and Morrison^) . The description of a streak 
lift-up (burst) as an "instability" initiated by the large-scale disturbances from the outer layer appears in 
the seminal papers of Kline et al££ and Kim et alJ^. MorrisonM^ describes more recent ideas concerning 
inner-outer interaction and its relationship to "inactive motion"^i2r^2i. 

Hall & Sherwir>i2£ take an alternative approach of considering inviscid waves in the wavy critical layer 
(where the wave speed is close to the convective velocity) of a streaky base flow. Using an earlier theory due 
to Hall & Smiths on vortex/wave interaction, they describe the nonlinear interaction of a self-sustaining 
process in which the nonlinear terms of finite-amplitude waves drive the streamwise vortices through a jump 
in the stresses at the critical layer, which in turn drive the streaks present in the base flow. The unstable 
equilibrium solution they determine describes the whole cycle of the self-sustaining process and can be 
loosely considered as a nonlinear eigenvalue problem. As such, it may be kept in mind where we use the 
description of a linear system driven by a nonlinear feedback forcing. 

Following Reynolds and Hussair>i2^, del Alamo and Jimenea^ have undertaken a temporal stability 
analysis of the Orr-Sommerfeld and Squire (OSS) equations in turbulent channel flow (i?e T = 2000) using 
a variable eddy viscosity: they show that maximum amplification of disturbances occurs at two spanwisc 
wavelengths, one corresponding to the widely accepted streak spacing, A+ 3 = 100, ( the '+' superscript 
denotes a variable non-dimensionalised by the viscous length scale, v ju r where v is the kinematic viscosity 
and u r , the friction velocity, u T — \J r w /p, where t w is the wall shear stress and p is the density), the other 
occurring at X X3 = 3h, where h is the channel half-height. While the former is clear evidence of near wall 
streaks (with streamwise wavelength of \ Xl « 1000), the latter indicates the presence of VLSMs. They also 
note that the fluctuations in streamwise velocity contain nearly all the kinetic energy and last longer than 
those in the wall-normal velocity. In a similar vein, Hwang and Cossu^ have shown that, in a turbulent 
channel flow for a sufficiently large Reynolds number, two distinct peaks of optimal growth appear, one 
scaling with viscous scales, the other with outer scales. 

McKeon and SharmaA^£ have explored a simple, essentially linear, forcing-response type description of the 
dominant processes in high- Reynolds- number turbulent pipe flow. The model reproduces inner scaling of 
the small scales close to the wall and outer scaling in the flow interior and displays features representative 
of VLSMs including their modulation of the smaller scale features. In contrast to Landahl's theory, the 
work has addressed scaling with Reynolds number. However, both theories stress the importance of linear 
mechanisms. One important feature of the theory is that a high response to forcing is observed around the 
critical layer, and in regions of high shear. At higher Reynolds numbers, the theory predicts that the effect 
on the near-wall turbulence of VLSM-type structures become more important, whereas at lower Reynolds 
numbers the lift- up mechanism is more important. This prediction is supported by the control simulations 
of Touber & Leschziner-^. 
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D. Passivity 

Passivity is an energy concept; its origin lies in circuit theoryi^ 9 -. A component is called passive if only a 
finite amount of energy may be extracted from it. To take a physical example, if f(t) is a forcing function 
or field at time t on a system Z and v(t) is its velocity, such that v = Zf, then the power consumed by 

the system from time t = to time T is (/, u)j T j = J Q f(t)v(t) dt. If we assume the initial conditions are 
zero, and this integral is positive for all T, then the system is passive. Essentially, it is a statement that 
the instantaneous power consumption is always positive. If Z is linear, the requirement for (/, v)r T i > is 
then equivalent to the requirement that Z is positive real, Z(jco) + Z*{juj) > 0. To show this, the integral 
{/, v), T i = (/, Zf), Q T j is considered in the frequency domain. 

The passivity theorem simply states that the feedback interconnection of two passive elements is itself pas- 
sive. Intuitively put, if two elements which cannot produce energy are connected in a feedback arrangement, 
then the feedback arrangement as a whole also cannot produce energy. We will show this for our particular 
case, and for general proofs and more information the reader is directed to standard control text o 79 ' 130 . 



E. Landahl's theory, scales and waves 



It is instructive to preface the current analysis with a review of the basic ideas of Landahl's theor y 102 ! 131 " -^, 
as inspired by the early wave theories of the viscous sublayer (Sternbergi 3 ^ and Morrison, Bullock and 
Kronauepi 3 ^) . Here we define a wave as a motion that exhibits a convection velocity that is constant over a 
region in wall-normal distance. An accepted definition of wave motion is one in which energy is transported 
but without bulk motion: hence the wave motion refers only to the fluctuating pressure and velocity field. 
Note that this is a stronger requirement (and a more physical definition) than a superposition of Fourier 
modes (Phillips 98 ). Therefore defining the viscous sublayer as a wave guide in which the least-damped waves 
exhibit significant correlation over large distances*^!, while useful, should not be taken too far. 

The current approach is reminiscent of Landahl's ideas, because both theories have the nonlinear terms 
as a right-hand side forcing to the linear problem as a common point of departure. However Landahl's 
theory considers perturbations to the mean turbulent profile, whereas we will consider perturbations about 
the laminar one. This difference will be examined more fully below. Specifically, Landahl's theor y 102 i 132 i 133 
considers the Orr-Sommerfeld and Squire equations, which respectively, may be written as 

DX7 2 u 2 -"du 2 V 4 u 2 



Dfj2 <u'^-^=r, (2) 



Dt dx 3 Re 
where 

D ( d - d \ 

D-t = {dt +U d^)> < 3 > 

U is the mean flow profile U(y)8n subject to a three-dimensional disturbance with velocity Ui(xj,t) and 
pressure, p(xi,t). Here the streamwise direction is x%, the wall- normal direction is x 2 and the spanwise 
directions is X3, giving the total velocity field as Ui(xj,t) = U(y)6u + Ui(xj,t). The wall-normal vorticity is 
i) = V x u. The forcing terms q and r are quadratic terms involving Reynolds stresses. Landahl proposed 
that q and r are significant only in localised regions in space and time, thus giving rise to "compact" source 
terms in (fTJ) and <j2j) (see for example Landahl 1 —). The picture of sublayer motion is therefore one in 
which regions of "intense small-scale turbulence of an intermittent nature" are interspersed by periods of 
"laminar-like but unsteady motion of larger scale" . 

Landahli 3 ^ identified three timescales associated with parallel mean shear flow, each a measure of the 
time after the creation of the structure from the original disturbance: 



G 




FIG. 1. Two-scale model of near- wall turbulence showing inner-outer interaction (after LandahU^). The small-scale 
wavelength is denoted by X s and the larger scale by A ( . 



1. the shear interaction timescale 



t s = U 



where t w — [iU w . Hence tf = 1; 
2. the viscous interaction timescale, 



tu = 



3. and the nonlinear timescale, 



t n =L/u , i+«100; 



where L and uq are the streamwise length scale and velocity scale, respectively, associated with the initial 
disturbance. Note that here t+ 3> tj, while conventional turbulence timescales require l/u <C I 2 /v where 
Re = ul/v is large. For short times after the creation of a structure from the original disturbance, the effects 
of both viscosity and nonlinearity may be neglected; in particular, the nonlinearity is assumed to operate 
only during short intermittent bursts of a local, secondary instability. Neglecting viscosity and linearising 
gives the Rayleigh equation for disturbances to parallel inviscid flow, 



DV 2 u 2 jj" du2 _ 

Dt dX! 



(4) 



In addition to the three timescales, Landahl described two scales of motion as important in understanding 
turbulent shear flow. Decomposing the velocity field into large- and small-scale components, 

Ui=<2? + fif, (5) 

the motions at the wavelengths of the large and small scales, A ( and A s respectively, are assumed not 
to interact if e = A s /A' « 1. A schematic of this description is given in Figure [TJ Substitution of ([5]) 
into Eqs.(UJ and @ and retaining only terms of leading order in e provides a pair of equations, one each 
for the large-scale field and the small-scale field of the form of ([TJ and ([2]). Later we discuss the relative 
importance of the different timescales and length scales to flow control and find these concepts useful even at 
relatively low turbulent Reynolds numbers. Landahl's scale separation has much in common with the ideas 
of Townsen d 119 ' 121 and BradshawiSIi concerning inactive motion and inner-outer interaction (Morrisonii^). 

This description also provides a basis for a more formal analysis in which wavenumber-frequency spectra 
are dominated by the least-damped Orr-Sommerfeld waves near resonance (see Landahl^i). For waves to 
be identifiable, it is necessary to form the wavenumber-frequency spectrum from which an unambiguous 
convection velocity can be obtained: then wave motion will appear as a reasonably narrow convective 
ridged. The response may further be localised at a wall-normal location and is stronger for certain mode 
shapeai^i. Barkis has identified the wave- like structure of the near- wall v,2— component motion as arising 
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from the Orr-Sommerfeld eigenvalues, while some eigenvalues relating to the the horizontal components 
correspond to the viscous decay of wall-parallel motion. Both Landahi^ and Russell and LandahU^ note 
that horizontal pressure gradients are small during these "quiescent" periods, and much smaller than those 
associated with a lift-up. 

The large- and small-scale decomposition raises the question of resonance. Jang, Benney and Gran-^i 
have proposed that, if the forcing function occurs at a frequency /wavenumber combination that matches 
the leading eigenmodes of the Squire equation, then resonant forcing occurs and there is the potential for 
large growth in amplitude before viscous damping occurs. They have shown that such a resonance could 
occur at a spanwise wavenumber corresponding to a streak spacing of 90 wall units. Zaki and DurbiiJ^ 
have discussed resonance in the context of the spatial problem and have shown that the dispersion relation 
for the homogeneous Squire operator is identical to that of the Orr-Sommerfeld equation, making resonance 
possible. However, Hultgren and GustavssoiJ^ have noted that since this growth mechanism is associated 
with the continuous spectrum, it is only possible when the flow is semi-bounded. Since our study is for a 
closed flow, we only encounter discrete modes. McKeon and Sharma offer an interpretation of resonance 
in terms of pseudospectrapi2£. Essentially, they understand this resonance as the high (but non-singular) 
system response to harmonic forcing resulting from left half-plane (stable) eigenvalues approaching the 
imaginary axis. Truly neutral or inviscid modes would be located at the imaginary axis. This high response 
is manifested as a high resolvent norm. 



II. MODEL FORMULATION 



This section describes the model formulation. We consider a three-component velocity field perturba- 
tion u(x,t) about an assumed time- independent solution Uq in the presence of a divergence- free, bounded 
exogenous disturbance forcing. This gives the net velocity vector field 

Ui = U 0i + Ui. (6) 

The steady pressure is similarly perturbed by p{x,t). 

In the next section, we will seek a control function f(x, t) to globally stabilise an assumed time- independent 
solution Uq in the presence of a divergence-free, bounded exogenous disturbance forcing d(x, t), representing 
unmodelled disturbances such as that arising from vibration, thermal disturbance, etc. This solution may 
or may not be stable to small perturbations in the absence of control. The turbulent mean is, in general, not 
a time-independent solution. Substitution of (|6|) into the Navier-Stokes equations gives the perturbation 
equations 

dUj / duj dllpi 
dt \ 3 dxj 3 dxj 
dp 

+ th — - 
i 

du 



h 3 



Uj dxi'- 



(7b) 



A substitution has been made for the nonlinear term, giving coupled linear and nonlinear equations. The 
pressure term is eliminated, along with the divergence equation, by the projection n onto the space of 
divergence-free functions. We do not make the linearising assumption of small perturbations. 

The forcing resulting from the control is restricted by a linear operator B(x), representing physical lim- 
itations on the actuation. The range of B spans the space of the divergence-free body forcing arising from 
all possible control actions. Thus, Bf(x,t) is the forcing on the fluid arising from the control function at 
time t and position x. For the purposes of understanding the current simulations, we may consider B as the 
identity when acting on the wall-normal velocity and zero otherwise. Let y(x, t) be the measurements made 
at time t, modelled by yi = CijUj, so that C{x) is a linear operator mapping the flow field to y. 
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A. The discretised equations 



Discretised, the equations ([7]) have the state-space form 

z =Az + Bin + B x d + B 2 f (8a) 
yi =Ciz (8b) 
y 2 =C 2 z. (8c) 

Matrices C\ and B\ are only used at the synthesis stage. They give respectively (C%) the flow field 
at the discretisation points weighted such that yi{t)'y\{t) approximates the perturbation energy E(t) = 
J x&n Ui(x,t)ui(x,t) dx, and (i?i) the forcing on the flow field from the nonlinearity, similarly weighted. To 
find C\ we require the mesh weighting appropriate for the discretisation chosen. For the discretised state 
the perturbation energy E(t) is approximated by the inner product on a positive-definite matrix R so that 
E(t) ~ z(t)* Rz(t), approaching equality in the continuous limit. Thus we require simply G\C\ — R. The 
input matrix B\ associated with the forcing from the nonlinearity n is determined similarly. The matrix C 2 
simply gives the measurements y 2 from the flow field state z, and is the discrete approximation of operator 
C. Similarly, B 2 is a matrix approximating operator B, describing the effect of the actuation on the flow 
field. 

State-space representations of the linearised Navier-Stokes equations are well known in the literature (see 
Bewley 85 for a pedagogical example). However, our formulation explicitly retains the nonlinearity as a 
forcing. 



III. CONTROLLER SPECIFICATION AND DESIGN 



This section specifies the requirements on the controller. For the following we will make use of the 
temporal-spatial inner product 

(a,b)= [ [ a*(x,t)bi(x,t) dx dt (9) 
Jo Jxen 

and the purely spatial inner product 

[a,6]= / a*(x,t)bi(x,t) dx. (10) 
The pair a and b is passive if (a, b) > 0. 

We define the perturbation energy as the spatial L 2 norm (induced by ()10|) ) of perturbations from the 
laminar profile, and the turbulence kinetic energy as the spatial L 2 norm of the perturbations from the 
turbulent mean profile. 

The rate of change of the total perturbation energy E of viscous shear flow integrated over a closed domain 
fl C M 3 is given by the Reynolds-Orr equation, which in parallel shear flow is 

f = -/ xen (W™ + ^) ^ (11) 

where U^y) is the y-derivative of the laminar profile Uo(y) and D is the dissipation rate. As before, x 
denotes a point in il. In essence, the aim of the controller is to provide actuation such that ^1 < for any 
disturbance. 

The broad design objectives of stability and robustness are achieved by application of the passivity theo- 
rem, which gives general, open-loop conditions for closed-loop stability of two arbitrary elements connected 
in a feedback loop. 
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FIG. 2. The feedback loop for controlled Navier- Stokes. The system inside the dashed box is Q. 



A. The nonlinear and pressure terms 



We will use the fact that the nonlinear term (|Tb[) is passive, specifically that (u, n) = for all T. Applying 
the divergence theorem and the divergence-free condition, it is easily shown that the inner integral in this 
expression is equivalent to an integral over the boundary dfl, 

UiUi dx = / (uj(x,t)uj(x,t))ui(x,t)£i dx (12) 

where £ is the outward-facing unit vector perpendicular to the boundary of the flow domain. 

Physically interpreted, f| 12|) quantifies the net flux of disturbance energy out of the domain through the 
boundary per unit time. In a closed or periodic domain, the contribution to this integral from volume 
forcing is necessarily zero. However, in an open domain, or with transpiration at the boundary, the flux of 
disturbance energy through the inlet and outlet boundaries and the net rate of flux of disturbance energy 
from any boundary control would both contribute. Were there such a contribution, (fT2|) would enter as 
a nonlinear constraint on the control law. For the open-domain case there will be a net flux out of the 
domain of the disturbance energy, where the outflow of perturbation energy is greater than in inflow of 
the perturbation energy (i.e. with relatively quiet inflow conditions). In these cases, the nonlinearity has a 
stabilising influence in the domain of study. In our case of periodic flow with body forcing, however, the 
integral is zero. 

Similarly, the contribution due to the pressure term is given by, 

Op f 

Ui— — dx — I pu^i dx — 0. (13) 
xen " x i Jxeon 

Again, in the case of actuation by boundary transpiration or in an open domain, this term would contribute 
to the perturbation energy. In our case, the integral is zero, simplifying the projection onto a divergence-free 
basis. 



B. The closed-loop linear terms 



In this section we will consider the discretised system equations (j5J) (the development is almost identical 
for the original perturbation equations (J7J). We can write the discretised equations in the compact matrix 
notation 
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where we have defined e = d + n. Let G be the linear system which, when discretised, has state-space 
realisation (|14[) taking e and / to y\ and y 2 . 

Furthermore, define K as the feedback law which, when discretised, generates the control action / from 
measurements y 2 , with state-space realisation 
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(15) 



Define Q as the system which maps the forcing from the nonlinearity to the fiowfield, u — Qe. Once 
discretised, Q is therefore the closed- loop of and (|15l) . This arrangement is depicted as a block diagram 
in Figure[2j with Q being the system inside the dashed box. Eliminating / and y 2 , the discretised Q therefore 
has a state-space realisation 
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(16) 



C. Stability 

The control problem therefore is to find a control forcing / such that u(x,t) — > as t —> 00, given 
the measurements and any exogenous bounded disturbance d. To do this, we must consider the stability 
properties of the system as a whole. Applying the passivity theorem, if N is passive and Q is strictly positive 
real, then the closed loop in Figure[3] (representing the controlled Navier-Stokes equations) is internally stable 
and is also strictly passive. In other words, (u,d) > 0. The case (u, d) = would imply that the forcing d 
acts orthogonally to u, and (u, d) > implies that d acts to reduce u in the feedback system. This result is 
simply verified; from Figure [2] and by the strict positivity of Q and the passivity of A/", 

(u,d) =(u,e-n) 

= (u, e) - (u, n) > 0. ' ' 

Note that if the uncontrolled, linearised plant is already passive, no control is required, as u is already 
bounded. The expression (u, d) quantifies the flow perturbation energy due to the disturbance. Physically, 
passivity of the controlled flow implies it only dissipates perturbation energy. Since this is true for all 
bounded d, it implies that all disturbances eventually decay. 

To achieve this, we wish to find a controller K such that the discretised Q is strictly positive real, or 
equivalently, Q(s) + Q*(s) > or (e,yi)j 0T j > for any T. This synthesis problem is a fairly standard 
exercise in robust control theory. The details of the synthesis procedure are in Appendix [X] Bounds on the 
perturbation energy production of the controlled flow are given in Appendix [Bj 



IV. APPLICATION TO PERIODIC CHANNEL FLOW 



To test the controller, we consider three-dimensional perturbations to plane Poiseuille flow at Re T = 100, 
with a constant mass flux. For this arrangement, laminar flow is the state with minimal sustainable drag, 
once control effort has been taken into consideratio n 91 ! 143 . Accordingly, we aim to sustain the laminar flow 
state at conditions where it would otherwise not persist. The flow domain is the space between two plates 
parallel in the (xi, X3)-plane, at x 2 = ±1. x\— and X3— direction periodicity is assumed. We consider 
actuation in the whole domain provided by simple body forcing of the wall-normal velocity which is also the 
measurement. 

The geometry allows Fourier transform of the linearised problem in the x\ and X3 directions which 
converts the spatially continuous problem into a number of decoupled continuous problems at particular 
Fourier wavenumber pairs. Truncation at suitably high wavenumber ensures an (x%, X3)-discrete problem 
with sufficient resolution. Further projection onto the Chebyshev polynomials in the wall-normal direction 
results in a number of linear time-invariant state-space control problems. The problem thus decouples at 
the (linear) synthesis stage, giving a block-diagonal A matrix in the state-space formulation. The decoupled 
wavenumbers only interact, via the nonlinearity, at the full simulation stage. The controller synthesis 
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FIG. 3. Pressure gradient variation with time: upper (-), uncontrolled; upper (- -) ki,k$ < 2; lower (- -), fci, &3 < 4; 
lower (-), fci, &3 < 8. 



problem is further simplified because control is unnecessary at the highest wavenumbers where we find that 
viscosity dominates and that the linearised system is already close to passive. 

Given that the source of the system non-normality is the interaction of the wall-normal velocity with the 
shear of the base profile, it is expected that actuation of wall-normal velocity only will be sufficient for the 
purposes of this study. 

The periodic spanwise and streamwise boundary conditions are naturally enforced by the Fourier trans- 
form. Further, any forcing is divergence-free as it is expressed in a divergence-free basis. 



A. Implementation 

The Reynolds number based on the target parabolic profile centre-line velocity and channel half-height was 
Re = 1709 (or Re T = 100 based on the friction velocity). The channel width was |7r and the length was An, 
which was deemed sufficient to provide accurate statistics based on previous studiesi^. The simulation was 
performed using a modified version of Channelflow 0.9. 15^^ which solves for the primitive variables using 
spectral discretisation in the spatial directions (Chebyshev in the wall- normal direction, Fourier otherwise). 
The flow field was advanced in time by a mixed third-order Runge-Kutta scheme, that treats the linear terms 
implicitly and the nonlinear terms explicitly. The number of modes used was 32 x 71 x 32. A variable time 
step was set capped at At = 0.01 which was sufficiently short to keep the CFL number low. The nonlinear 
terms were computed in skew-symmetric format with 3/2 dealiasing in the wall-parallel directions. The 
statistics of the unmanipulatcd flow were verified by comparing them to a database provided by Kuroda 
and Kasagii^. The profiles of the mean velocity and the Reynolds stresses largely collapse, with a slight 
discrepancy in the Mi-component normal stress near the centre-line. Given the accuracy, the discretisation 
is considered sufficient for this study. 

The control action was integrated using a zero-order hold, stepped at At = 0.01. The control penalty was 
set at £ = 0.01 as defined in (|A3[) . The control action was restricted to forcing on the wall-normal velocity 
and the sensing was likewise restricted. A value of 7 < 1.02 was achieved for the Cayley-transformed system 
at all wavenumbers, indicating that the controlled flow is very close to passive. 

Figure [3] shows the pressure gradient for the controlled cases decreasing with time. Interestingly, the 
restriction of control to lower wavenumbers (fci, £3 < 4) produces forcing that is almost as effective as that 
for the larger wavenumber range (fci , &3 < 8). The control fails when restricted to fci, < 2. This suggests 
that the effect of viscous damping is significant enough at the highest wavenumbers to overcome the energy 
production due to the shear interaction. Consequently, we infer that the dominant production mechanisms 
occur at these larger scales: the key requirement is that the control scheme should resolve streaks and 
streamwise vortices. Hence, additional control at ki,k 3 > 8 is ineffective because the forcing appears at 
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FIG. 4. Mean-square forcing (averaged over wall-parallel planes) at t — 10 (-) and t = 50 (• • ■). The forcing decreases 
substantially over time. It peaks close to X2 + ~ 20. 
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(c) forcing at X2 = h 

FIG. 5. Contours at various wall-normal distances of the forcing provided by the controller at t — 10 

scales shorter than the streak spacing (k\ ~ 8 corresponds to A3 + ~ 75). 

Figure 2] shows the mean-square averages over wall-parallel planes of the forcing as it varies with wall- 
normal distance at various times. Where variables are expressed in viscous units {e.g. x% + ), the relevant 
viscous scale is calculated from the uncontrolled flow (otherwise it would change with time) . The forcing is 
concentrated around X2 + — 20 where the shear interaction is most significant. The forcing peak decreases 
over time and moves further into the flow interior. This indicates that, as the laminar profile is approached, 
only minimal control effort is required. Contours of the forcing at various wall-normal distances are shown 
in Figure O 
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FIG. 6. Evolution of the average of (-), u| (— •), w| (■ ■ •) and p 2 (- -) with control at k\,kz < 4. The log scale 
shows that the rate of decline is fastest for U2 and p. 
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(a) contours of u\ at X2 + = 20, t = 10 
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(b) contours of u\ at = 20, t = 20 
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(c) contours of ui at X2 + = 20, t = 50 
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FIG. 7. Contours of ui for the controlled flow at a; 2 + =20, t = 10, 20, 50 



Next we examine the results in relation to the pressure field. Figure |6] shows the energy of the velocity 
and pressure fields over time, for one controlled case. Figures [71 HI HI and [TU] show contours of these fields 
at X2 + = 20 for various times. Comparison shows that the wall-normal and pressure perturbations are 
controlled very quickly. The spanwise perturbations subsequently decay, then lastly the energetic streaky 
streamwise contours. There is a brief spike in the pressure field as the controller comes on-line. This ordering 
supports the picture that the interaction of wall-normal motion and shearing is a minimum requirement for 
the control of streaks—. 
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FIG. 9. Contours of u 3 for the controlled flow at x 2 + = 20, t = 10, 20, 50 
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(a) contours of p at x^^ = 20, t = 10 
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(b) contours of p at x^^ = 20, t = 20 




(c) contours of p at = 20, i = 50 
FIG. 10. Contours of p for the controlled flow at x 2 + = 20, t = 10, 20, 50 

V. THE ROLE OF PRESSURE AND WALL-NORMAL VELOCITY FLUCTUATIONS 



The Poisson equation for pressure fluctuations in reduced form appropriate for channel flow is given by 



0X1 



c) 2 



dx\dx2 



U1U2 - U\U2 



(18) 



where we are once again considering perturbations to the turbulent mean profile. The first term on the 
right-hand side is the linear or "rapid" source and the second term is the nonlinear or "slow" term, the 
physical distinction coming from the fact that the linear term changes as soon as the mean rate of strain 
changes. Pressure fluctuations are known to be well correlated across a shear flow: Kim 1 — has shown that 
the two-point correlation extends over both large wall-normal distances and over large spanwise distances 
near the centre-line of turbulent channel flow. He has also shown that contributions come mainly from the 
slow source term, except close to the wall where contributions from the rapid and slow terms are about the 
same. The earlier measurements of Sternberg^ also indicate that the linear pressure fluctuation field at 
the edge of the sublayer is larger than the nonlinear field. 

The rapid term arises from the inviscid, linear term, associated with the interaction of the wall-normal 
disturbances with the mean profile. The proposition that this is the leading term that must be controlled, 
is consistent with the controller's success, and the localisation of the control near x 2 + — 20. This accepted, 
the control may reasonably be expected to work when restricted to actuation and sensing on wall-normal 
velocity or pressure disturbances. Landahl 1 ^ has formulated the perturbation field in terms of the pressure 
as: 
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(19) 



where s is the nonlinear term corresponding to that in equation (JlJ . Interpretation of the the terms on the 
left-hand side suggests that terms like V 2 p and V 4 p will be significant where p changes rapidly. 

More physically, Kim 1 ^ shows that, in channel flow at Re T = 179, the static pressure is only slightly 
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negatively skewed, but has flatness factors that are typically twice the Gaussian value of three over much of 
the channel height. Writing the mean-square acceleration as 



Dt J \dxj \dx] 




(20) 



Batchelor and TownsencU^ have shown that, at Reynolds numbers high enough for local isotropy (such that 
the diffusion term is negligible), the mean-square pressure gradient is much larger than the mean-square 
viscous force. Further, they suggest that 



20- 2 (^ , (2D 



where the constant is determined empirically. Dunn and Morrisoni^ - (see Figure ITT]) show that, outside the 
viscous sublayer, the factor is about 5 — 10, even at low Reynolds numbers. In the current work (Re T — 100), 
we observe it to be about 2 at x 2 + = 20 for the uncontrolled case, as shown in Figure IT2"1 

Equation (f20| suggests that the mean-square acceleration comprises prolonged viscous intervals 'pulsed' 
periodically by the mean-square pressure gradient, as illustrated in Figure [T51 Thus a pressure field distri- 
bution of small skewness, but large flatness, gives rise to a pressure-gradient distribution of which the first 
moment is very small, but with even moments that are significantly larger. Kimi^Z also shows that contri- 
butions to the mean-square wall pressure are principally local in nature even though the instantaneous wall 
pressure receives significant contributions from the opposite wall of the channel. Therefore, the mean-square 
pressure close to the surface is intimately related to the structure there. In terms of the sublayer populated 
with quasi-streamwise vortices, this means merely that a low-pressure region (approximately coinciding with 
the vortex core) always has two opposite-signed pressure gradients in the cross-sectional plane of the vortex. 
Hence, at x 2 + ~ 25, {dp/dxi) 2 rs 0.b(dp/dx 2 ) 2 ~ 0.5(dp/dx3) 2 . Kimi^l also notes that dp/dx is not a good 
indicator of quasi-streamwise vortices whereas the vertical and spanwise gradients are. 

The relevance of Landahl's equations lies in the fact that, over the short time for which the controller is 
active, the considerably longer turbulence timescale means that the turbulence itself is not very significant. 
They therefore offer an explanation of the controller's success. The controller reduces the pressure gradients 
over time through action on the v— component and the linear source term in equation (|18l) . In the short 
term (t < 35), the pressure term is higher than in the uncontrolled case, reflecting the controller's action. 
In the very short term, there is a brief spike, perhaps due to the initialisation of the controller. The viscous 
term decreases almost monotonically (not shown). The net result is that the ratio of the pressure terms to 
the viscous terms increases for a short time, then declines as the controller action takes effect. Figure [TTJ 
shows that the location of the maximum mean-square pressure gradient initially occurs at x 2 + ~ 20, where 
the forcing is a maximum. The controller moves the peak location of this term over time, presumably as 
the effective Reynolds number drops. 



VI. CONCLUSIONS 



A new characterisation of stabilising feedback laws for incompressible Navier-Stokes flows has been pre- 
sented in terms of passivity theory. The control is designed to make the magnitude of any perturbation to the 
laminar operating point decay monotonically. The flow equations are discretised and the ensuing controller 
synthesis problem results in two game-theoretic algebraic Riccati equations. When these Riccati equations 
have solutions, a globally stabilising, linear, controller can be synthesised. A simple synthesis procedure has 
been presented which is derived by the Cayley transformation of the positivity problem into an auxiliary 
7-optimisation T-Loo problem. Tools for the 7-optimisation problem are widely available in packages such as 
MATLAB®J£i or Octave^. A control effort penalty and measurement noise model has been introduced 
to avoid a singular control problem and its associated large control signals. The methodology allows an 
attempt at control with limited or insufficient actuation or sensing and permits bounds on the maximum 
perturbation energy production. It has been applied to turbulent channel flow with wall-normal interior 
body sensing and forcing. It was verified that the control relaminarised the flow, even when the forcing was 
confined to low wavenumbers (k x ,k z < 4). It seems likely that an important requirement is for the mean 
streak spacing to be resolved. Intuitively, this requirement explains the targeting of the shear interaction 
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FIG. 11. Mean-square pressure gradient and mean-square viscous force, defined by Equation (I21[l . but with the lower 
two figures using the target laminar profile, not the mean profile. The peak of the pressure term occurs at X2 + — 20 
until the control acts. 
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FIG. 12. Evolution of the ratio of the mean-square pressure perturbation gradient to the mean-square perturbed 
viscous terms at X2 + = 20 (uncontrolled (-), controlled ( — •) at k x ,k z < 4). 
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FIG. 13. Schematic of variation of mean-square acceleration with time: prolonged viscous periods pulsed by pressure 
'spikes' generated principally by quasi-streamwise vortices. 



mechanism. 

We have seen that controlling the wall-normal perturbations successfully, and with it the pressure pertur- 
bations, resulted in the eventual collapse of the streamwise streaky structures. This causality shows that 
the interaction of wall-normal motion and shearing is necessary for the formation of these streaks. 

The success of the control may be understood in terms of Batchelor and Townsend's result showing the 
importance of pressure-gradient fluctuations, and several essential features of Landahl's model. We observe 
that the shear interaction timescale is shorter than the viscous and nonlinear (turbulent) timescales. Since 
the shear interaction process is essentially linear and underpins the turbulent fluctuations, our control strat- 
egy is also linear. The shear interaction is governed by the wall-normal disturbance, which is related to the 
pressure via the linear ("fast") source term in the Poisson equation for pressure fluctuations. Consequently, 
the control may be satisfactorily restricted to wall-normal velocity or pressure. 

The response of the pressure and wall-normal velocity disturbances is particularly high for large, wave-like 
motions at close to the convective velocity, which correlate over significant distances in planes parallel to the 
wall. These waves are remarkably non-dispersive, with an approximately constant phase velocity. At higher 
Reynolds numbers, we might expect the effect on the near-wall turbulence of superstructures to become 
more important, with correspondingly more stringent requirements on actuation authority. These are also 
wave-like disturbances, but associated with a near-singular response of the linear terms to the nonlinear 
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FIG. 14. A flow diagram giving an overview of the entire controller synthesis procedure. 



convection termini which may be understood in relation to critical layer theory. 

In Landahl's theory, when the phase velocity is equal to the group velocity of disturbances emanating from 
further upstream, a secondary instability in the form of a burst occur o 133 ' 153 . Inhibiting the propagation 
of these waves precludes the occurrence of such nonlinear secondary instabilities associated with turbulent 
flow. Again, this aspect to the control problem is essentially linear and inviscid. 

Thus we may interpret the control via Landahl's theory; the linear mechanisms on the "rapid" timescale are 
most important, which is why the controller concentrates on the shear instability. The controller concentrates 
on regions where this mechanism is most active, acting to "normalise" the response of the system. This can 
be achieved through manipulation of just the wall-normal component of the velocity field. 
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Appendix A: The passivity control synthesis procedure 
1. Overview 

This appendix details the process of finding a discrete controller that satisfies the closed-loop passivity 
requirement for the discretised system. 

The aim is simply to provide the reader with enough information to replicate the procedure. No proofs are 
presented, and lengthy explanation is avoided. Such proofs are available in the book by Green & Limebeer— , 
which our presentation follows closely. That text, alongside Kailathi^, serves as a good background reference 
on the systems theory and robust control theory used in this work. 

This appendix is split up into a number of sections. The first section gives an overview of the synthesis 
procedure. The second section presents the transformation of the positive real synthesis problem into a 
general T-Loo problem. Then, the third section gives the loop-shifting transformations needed to convert 
the general problem into a simplified problem. Finally, in the fourth section, the solution to the simplified 
Tica control problem is given. 



2. Overview 



There are many possible approaches to solving the passivity control problem. In the approach chosen 
for this work, the synthesis problem is solved by applying a transformation to the system we wish to make 
passive, giving a new control problem where we have to find a controller to make the infinity norm of closed- 
loop with the transformed plant less than 1—. The resulting general problem is in turn solved using 
loop-shifting transformations^*^ and Riccati-based state space method o 79 ' 80 . The multi-step process is 
outlined in Figure [H] 
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The chosen approach is not necessarily the simplest (for instance, that of Sur>i££ is more direct), however 
it is robust and enables utilisation of readily available software such as the MATLAB® robust control 
toolbo»i^i. 

Four assumptions are made. Of these assumptions, one requires the stabilisability and detectability of 
the uncontrolled system. This is automatically satisfied if the flow is below the critical Reynolds number, 
where the first unstable eigenvalue appears. A second assumption is imposed to prevent unbounded control 
signals. The remaining assumptions are required for the solution method of the Riccati equations, and may 
be relaxed. 

The iterative method presented is useful in the case where there is insufficient actuation or sensing to 
make the closed-loop passive. In this case, we iteratively search for a controller to get the system as close to 
passive as is possible. This relaxation loses the strict guarantee of nonlinear stability, but it is still possible 
to quantifiably limit the perturbation energy production. We choose this method for our study, because 
it is more amenable to such a relaxation and so may be more applicable to cases with physical or design 
constraints on the available measurement and actuation. 

The solution to the general %oo problem is somewhat intricate, however it can be simplified using loop- 
shifting transformations summarised in Section [A 41 so that the simplified 1-Loq theory presented in Section 
IA 51 can be applied. This synthesis method requires the solution of two Algebraic Riccati Equations ( AREs) 
at each wavenumber. 

For brevity, the notation in each appendix is self-contained. 



3. Transformation of the positive real synthesis problem to a general T-Loo synthesis problem 



Let G have state-space matrices given by 
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The closed-loop transfer function of G and a controller K will be strictly positive real if and only if the 
closed- loop transfer function of G and K has infinity norm less than 1 (see Safonov et alM-), where 
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The problem therefore has become to find a controller K to minimise the norm, 7, of the closed-loop 
of K with G. If 7 < 1 , the closed- loop between the original system G and K is strictly positive real, thereby 
solving the original passivity problem. There is no a priori way to find a minimal 7, so it is necessary to 
perform an iterative search over 7. 



a. Control penalty 



The control problem as presented above permits unbounded control signals, essentially because it does not 
penalise the control effort. This is tackled by introducing a penalty on the control, and a model for sensor 
noise. The penalties are made orthogonal to the dynamics, by augmenting G (to give G + ), with scalar the 
penalty weighting e, 
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It will be seen that the penalty is necessary to satisfy rank assumptions on D\ 2 and D21 of Section 
At this point, we apply the results of sections I A 41 and lA 51 to find the required controller. 
We find that the control penalty is the primary obstacle to minimising 7. 



(A3) 
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4. The loop shifting transformations 



This section describes, without explanation, the loop shifting transformations required to convert the 
problem into a form that is solved for in Section lA 51 
From the previous section, we have a system of the form 





' A 


Bi 


B 2 


P(s) = 


Ci 


D u 






c 2 


D 21 


D 22 



(A4) 



with A G C nxn , B\ G C" xm , C 2 G C qxn , and other matrices dimensioned accordingly. 

The Hoc controller synthesis formulae are greatly simplified by assuming that Du — and D22 — and 
that D12 and D 2 \ satisfy some simple rank assumptions. 

Let the state-space matrices satisfy the following assumptions, 



Al. (A, £$2,02) is stabilisable and detectable, 
A2. rank(Z?i2) = m and rank(Z?2i) = q 
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The aim is to replace the system P with an equivalent problem involving P, where 
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with the simplified assumptions, 

Al. (A,B2,C2) is stabilisable and detectable, 

A2. D* 12 Di2 = I m and D 2 iD* 21 = I q 
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a. Minimise \\Dn\ 
We define 
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Begin by choosing F such that ||-Dn|| = 70 is minimised, where 70 = max{||i)j 2 £)ii||, H-Dii-DJill}- This 
can be done in more that one way (see Green & Limebeer— for further details) . 
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b. Eliminate Du 

Define 
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We can eliminate Du, by substitution we see directly that 
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c. Eliminate D22 

Eliminate _D 22 by connecting —D 2 2 in parallel with P 22 . 



d. Rank conditions on D12 and D21 

Find scaling matrices Si and S 2 such that Di 2 = Di 2 S\ with Dl 2 D\ 2 = I m , and similarly D 2 i = S 2 -D 2 i 
with Ai-DIi = /g. The rescaled system is then 
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e. Controller synthesis 

Find the controller K to solve the small gain problem for the system P in (|A9I) , using the method presented 
in Appendix I A 51 



f. Reversing the loop shifting 

The final step is to apply the preceeding steps of this appendix to the controller K in reverse. Where 



K = 





Bk 


_ Ck 






(A10) 



the final controller is given by 



K 



' A k + BuSxb^I - 


f FD22) 


1 S2C'k 


BkS\ — B k SiD 2 2{I 


-Db h b 22 )- l F~ 


S 2 C k -FD 2 2(Il 


-FD 22 y 


S 2 Ck 


-F + FL> 22 (/4 


-FD 2 2)- l F 



(All) 
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5. Solution to the simplified Hoo control synthesis problem 



In this section we present the controller synthesis formulae solving the small gain problem of Section [A4el 
above. 

Suppose the system P, given by 





' A 


B 1 


B-2 ' 


P(s) = 


Ci 





P12 




c 2 


D 21 






(A12) 



satisfies the simplified assumptions of Section \A 41 We seek a controller K such that the closed-loop of P 
and K is stable and the infinity norm of the closed- loop is less than 7. 
There exists such a K if and only if 



1. There exists a solution X to the ARE (|A13|) such that A - {B 2 B 2 - j^BxB'^X is asymptotically 
stable and X > 0. 



2. There exists a solution Y to the ARE (|A14j) such that A - Y{C 2 C' 2 - 7 _2 CiC() is asymptotically 
stable and Y > 0. 

3. The spectral radius, p(XY) < -f 2 . 
The AREs in question are 

XA + A'X - X{B 2 B' 2 - ^BxB'^X + C[(I- D 12 D' 12 )d = (A13) 
with A = A - B 2 D' 12 Ci, and 



AY + YA! - Y(C' 2 C 2 - 7- 2 C(Ci)Y + B X {I - D' 21 D 21 )B[ = 

with A = A - BiD' 21 C 2 . 

When these conditions are met, one such controller is given by 



K 



\A k 


B k ' 


[c k 






(A14) 



(A15) 



with 



A k = A + 7 - 2 Bi J B 1 A - B 2 D' 12 d + B' 2 X 

+ B k (C 2 + 1 - 1 D 21 B[X), 
B k = B 1 D' 2l Y{I - 7 - 2 Ay)- 1 (C 2 + ^D^X)', 
C k = —D' 12 Ci — B 2 X. 



Appendix B: Bounds on the perturbation energy 

The solution of the auxiliary small-gain problem (7 < 1) results in monotonic decay of the disturbance 
energy. In the case that 7 > 1, this property may be lost, however the method does optimise for the 
worst-case perturbation energy production. This is seen from the following argument. 

A transfer function Q is strictly positive real, if and only if its Cayley transform Q has infinity norm less 
than 1, i.e. ||Q||oo < l 95 - The Cayley transform Q of system Q is given by 



Q( s ) = (Q( s )-I)(Q( s ) + I)-\ 



(Bl) 



We have transformed the problem of making some transfer function Q(s) as close as possible to positive 
real into an equivalent problem of making Q(s) bounded real, i.e. ||Q||oo < 7- Then 



det[7 - 7 _1 Q(s)] ^ 0, for Re(s) > 0. 



(B2) 
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Using the Cayley transform (|B1[) it is straightforward to show that 



Q(s)Q*(s) 

=(Q(s) - I)(Q(s) + iy'iQ^s) + - I) 

<7 2 7. (B3) 



Rearrangement gives 



()(,) - : ( .s) > l—-L(Q*(s)Q( S )+I). (B4) 



As 7 — > 1, Q(s) becomes positive real. Bounding the right hand side by —a, 

■1-7 2 



-a = inf 



i + r 



(Q*(s)Q(s) + I) 



means a > (since 7 > 1). 

If u = Qe, then it is straightforward to show 

(u,e)>~(e,e) Ve. (B5) 

Since ||w||2 is the perturbation energy, this bounds the rate of perturbation energy production by any 
disturbance e and optimising 7 optimises this bound. 
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